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Abstract 

Recent developments on the study of mixed morphology supernova 
remnants (MMSNRs) have revealed the presence of metal rich X-ray emit- 
ting plasma inside a fraction of these remnant, a feature not properly ad- 
dressed by traditional models for these objects. Radial profiles of thermo- 
dynamical and chemical parameters are needed for a fruitful comparison 
of data and model of MMSNRs, but these are available only in a few 
cases. We analyze XMM-Newton data of two MMSNRs, namely IC443 
and G166.0+4.3, previously known to have solar metal abundances, and 
we perform spatially resolved spectral analysis of the X-ray emission. We 
detected enhanced abundances of Ne, Mg and Si in the hard X-ray bright 
peak in the north of IC443, and of S in the outer regions of G166.0+4.3. 
The metal abundances are not distributed uniformly in both remnants. 
The evaporating clouds model and the radiative SNR model fail to repro- 
duce consistently all the observational results. We suggest that further 
deep X-ray observations of MMSNRs may reveal more metal rich objects. 
More detailed models which include ISM-ejecta mixing are needed to ex- 
plain the nature of this growing subclass of MMSNRs. 



1 Introduction 

Mixed-morphology supernova remnants (MMSNRs) have been tradition- 
ally defined as remnant having a shell morphology in radio and a cen- 
trally pea ked morphology in t he X-ray band characterized by a thermal 
spectrum l|Rho fc P ctrc 1998). The origin of this atypical morphology 



is controversial. Tra ditional models mostly relies on the effects of ther- 
mal c onduction (e.g. ICox et al.lll999l . IShelton et al.lll999l , IWhite fc Lond 



Il99ll ). but there are also models which invoke projection effects ( Petrukl 
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120011 ). iTillev et al.l l|2006T ) tried to refine previous analytical models by 
using a set of hydrodynamical simulations aimed at reproducing the mor- 
phology observed in MMSNRs. 

In parallel to these theoretical efforts, there is a growing observational 
interest around MMS NRs, which has led in a few cases to useful com- 
parison with models. ISlane et"afl l|2002h favor the evaporating clouds 
IWhite fc Lond dl99ll) mode l for G290. 1-0.8 (rather then the radiative 
model of IShelton et al.lll999l ). b ut problems remain in th e comparison of 
the surface brightness profiles. lLazendic fc Slanel (|2006h also find a rea- 
sonable agreement between both models and X-ray observations of CTB1 
and HB21, but they point out that the predi cted central density are b elow 
the derived densities of the emiting plasma . iRho fc Borkow ski (2002]), in- 



stead, argue that the I White fc Long (| 1 99 lh evaporating cloud model is at 
odds with the presence of a strong radiative shock front in the MMSNR 
W28, but also the radiative model is at odd with th e strong thermal varia- 
tion found in this remnant (see also lChe valicr 1999 for a critical approach 
to the evaporating cloud model of MMSNRs) . 

In addition, several authors have pointed out that part of the X-ray 
emission of MMSNRs may be composed by thermal plasma with high 
metal abundances, as those observed i n stellar ejecta ins ide young his- 
torical SNRs (e.g. IShelton et all 12004 IChen et alj l2008h. or in other 
Gala ctic (e.g. Vel a SNR. iMiceli et al.l 20081 ; Cygnus Loop. lKatsuda et all 
120081 : Pup pis A. Illwang et al. l l200Sf> and Magellanic Cloud SNRs (e.g. 
0103-72.6, Park et al.ll2003al ; N49B. iPark et al.ll2003bl'). This im portant 
result was addressed and summarized by Lazendic fc Slanel (|2006T l ■ which 
report a compilation of 26 MMSNRs, 10 of which show some sign of en- 
hanced metal abundances in the X-ray spectrum. It is important to stress 
at this point that the presence of enhanced metal abundances inside MM- 
SNRs has not been addressed in detail yet in the models developed for 
this sub-class of remnants. In general, published models with thermal 
conduction do not take into account the mixing of ejecta with the cir- 
cumstellar and interstellar medium, so the emerging class of metal rich 
MMSNRs is still not properly understood. This is also the reason why 
the comparison between models and observation has mostly focused only 
on morphological issues, with the general comment that the presence of 
the additional ejecta component may reconcile the discrepancy betwee n 
model and observed profile (like in the case of W44. IShelton et alj|2004h . 

The mixed morphology category of SNRs (and especially the emerg- 
ing subclass of metal-rich MMSNRs) can be better understood if we have 
a large sample of objects with well studied characteristics. In particular, 
the accurate determination of density, temperature and abundances profile 
for MMSNRs will allow us to perform an accurate comparison with both 
traditional models and a new magnetohydrodynamical (MHD) model of 
shocked ISM and ejecta inside MMSNRs, subject to a forthcoming publi- 
cation. 

We have unde r taken a systematic review of the MMSNRs listed in 
lLazendic fc Slanel | 2006) with the primary goal to measure the abun- 
dances in the inner region of the remnants and other X-ray properties. 
In this paper, we review the XMM-Newton archive observations of two 
bright MMSNRs, IC443 and G166.0+4.3, with the aim of deriving the 
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distribution of temperature and metal abund ances inside these r emna nts. 
Both our targets are mentioned in the list of lLazendic fc Slanc (20 0j|) a s 
M MSNRs with s tanda rd metal abundances. However, Troia et al.l (2006) 
and lTroia et al.l (|2008T ) have studied in detail the X-ray emission of IC443, 
and they have found enhanced metal abundances in a remarkable ring fea- 
ture around the pulsar wind nebula, in the southern part of the remnant, 
attributed to emission from shocked ejecta. The compilation of MMSNRs 
with enhanced metal abundances needs therefore to be updated with new 
results from deep X-ray observations of MMSNRs. On the light of these 
new observational results, traditional models should be further verified 
and revised as necessary. 

The paper is organized as follows: in Sect. 2 and 3 we present the 
results obtained form XMM-Newton archive observations of IC443 and 
G166.0+4.3, respectively. In Sect. 4, we present the comparison of our 
results wit h the standard mode ls for MMSNRs, namely the evaporat- 
ing clouds IWhite fc Lond (Il99ll ) model and the radiative SNR model of 
ICox et al.l(l999l ). In Sect. 5, we present the conclusion of our work. 



2 IC443 

IC443 has been observed in several bands and extensively studied in the 
literature and, for this reason, it is a case study for the interaction be- 
tween SNR shock s and molecular clou d s. Recent X - ray st udies includes 
iTroia et al l (120081 ), iBvkov et all ||200S|). iTroia et al l (l2006h . while for in- 



frared and radio studie s see Noriega-Crespo et alj (12006 ^_ Bvkov et al.l 



20081 ), iLee et all j|2008l ). iRosado et al.l (|2007l ). iLeahvl (|2004D and refer- 



en ces therein. The re mnant has been classified as mixed-morphology 
by iRho fc Petrel (Il998l), based on observations with the ROSAT satel- 
lite. ITroia et all i|2006T ) and ITroia et all (|2008l ) confirm the centrally 
peaked morphology between 0.5 and 5 keV, with a bright X-ray peak near 
06 h 17 m 09 3 +22°45 m in the 1.4-5 keV (Fig. 0. However, they note that, 
at very soft energies (0.3-0.5 keV), the remnant shows an incomplete shell 
extending from N to E, corresponding to one of the region of interaction 
with a large cloud. They have also investigated the spatial distribution 
of the metals inside this remnant, using the technique of the equivalent 
width images. According to them, the distribution of Si and S is far from 
being homogeneous. In particular, a ring of metal rich plasma surround- 
ing the pulsar wind nebula is evident in the equivalent width images of 
these two elements (see their Fig. 3). Moreover, the northern region of 
the remnant, coincident with the maximum of the hard X-ray image (Fig. 
[!}, is also characterized by a large value of the Si and S equivalent width, 
suggesting higher metal abundances with respect to surrounding regions. 
In the following, we focus on this region. 

The data we use in this work are part Calib ration and Perform ance 
Verification phase of the XMM-Newton satellite l|jansen e t al. 200l|), for 
a total of 6 observations of IC443. In addition, we have used the public 
archive observat ion 0301960101. T he datasets we have used are the same 
datasets used bv lTroia et al.l (120081). We screened both PN and MOS data 
using the algorithm suggested bv lSnowden fc Kuntj l|2007l) . 
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Figure 1: Top: XMM-Newton EPC mosaic of the northern region of the IC443 
supernova remnant in the 1.4-5.0 keV band. The image is background sub- 
tracted and vignetting corrected. The pixel size is 10 arcsec. Five contours 
levels at 8, 10, 12, 14 and 16 x 10~ 4 cnt s -1 pix 1 are overlaid. These contours 
have been used as boundaries for the 5 spectral extraction regions discussed in 
the text. Bottom: same region at 1420 MHz (total intensity map) with X-ray 
contours overlaid (adapted from lLeahy||2004l ). 
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The results are reported in Fig. [3] In order to study the profile of the 
thermodynamical and chemical parameters of this remnant, we have cho- 
sen 5 spectral regions denned in terms of contours of iso-surface brightness, 
centered in the bright northern region, as shown in Fig. [T] For each region, 
we have computed a mean distance from the X-ray peak as the average of 
the distance of all the pixel in that region. XMM-Newton PN and MOS 
spectra have been extracted using the software SAS 7.1 and co r rected 
for detector non-uniformity using the task evigweight. iTroia et all (|2006l ) 
have shown that the X-ray emission in the IC443 interior is composed 
by two thermal components, the soft one with temperature in the range 
0.3-0.7 keV, in non-equilibrium of Ionization (NEI), and the second one 
with kT = 1.1 — 1.8 keV, in equilibrium of ionization (CIE). The soft com- 
ponent has been associated with the shocked interstellar material, while 
the hot component to the shocked stellar ejecta. Theref ore, we have 
used the same emissio n models, namely the mekal model l|Mewe et alj 
Il985ll and the VNEI dBorkowski et all l200ll ) model of XSPEC vll.3.2 
(|Arnaudlfl996h . We have fixed the chemical abundances of the soft com- 
ponent to solar values, and left the abundances of the hot component 
free to vary. The adopted models represent an approximation of the true 
conditions inside this complex remnant. However, this approach allows 
us to measure quantitatively the emission measure weighted temperature 
and metal abundances for each component in each spectrum. Moreover, 
by using the region layout shown in Fig. [T] we may derive meaningful 
constrain to any trend of the measured quantities related to the radial 
distance of the regions. 

In Fig. O we show an example of spectrum, while the complete results 
of the spectral analysis are reported in Fig. [3] for the 5 spectral regions, 
displayed in order of increasing distance from the putative center. The 
temperature does not show large variations, indicating that thermal con- 
duction must be efficient. The innermost regions are characterized by 
a higher metal abundances for Ne, Mg and Si. In these cases, the pro- 
files show a decreasing trend with increasing distance form the center. 
While Ne and Mg retain over-solar of abundances in most regions, Si is 
above 1 only in the two innermost regions. On the other hand, S and Fe 
seems to have more uniform abundances, about 1.0 and 0.3 respectively. 
The thermodynamical parame ters of both the co mponents are in general 
agreement with the findings of lTroia et alj (|2006T ) . T he abundances pat - 
terns are generally consistent with the ones found by ITroia et all (2008) 
in the south region of the remnant, but smaller in value if compared 
to the re markable ejec t a ring surrounding the pulsar wind bebula found 
by them. ITroia et alj |2008) have also compared the Mg/Si, S/Si and 
Fe/Si ratios of the ring w ith the nucle osynthesis model of core-collapse 
(jWooslev fe Weaver]! 19951 ) and Type la (jBadenes et aUl2003h supernovae, 
finding more agreement with the former. We have verified that Mg/Si, 
S/Si and Fe/Si ratios found by us in the innermost spectral region in Fig. 
[T]are in agreement with those of the ring region, thus we concluded that 
our results also support a Type II progenitor for IC443. 
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IC443 Third contour region 
PN (black), MOS1 (red) and MOS2 (green) 
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Figure 2: XMM-Newton EPIC spectrum of the third contour region in Fig. [T] 
(average distance 2.8' in Fig. [3}. The best fit two thermal component model 
is shown as continuous line (with residual in lower panel), while individual 
components are shown in dashed black lines for the PN spectrum only. Emission 
lines from O (0.6-0.8 keV), Ne (~ 1 keV), Mg (~ 1.4 keV), Si (1.8-2 keV), S 
(~ 2.5 keV) and Ar (~ 3 keV) are visible. 
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Figure 3: Results of the spectral fittings to the hard thermal component ob- 
tained in the 5 iso-brightness regions in IC443. See Fig. [I] for the position of the 
contours levels used for the region definition. Abundances are relative to solar 
values. The abscissa contains the average distance from the center (horizontal 
error bars are standard deviations). Vertical error bars are 2a uncertainties. 
The normalization of the spectrum is expressed in terms of n 2 l, where n is the 
plasma density and I is its extension along the line of sight. 



3 G166.0+4.3 (VRO 42.05.01) 



G166.0+4.3 is a supernova remnant in the anti-center direction whose 
X-ray and radio morphology has d rawn considerabl e inter est in the liter- 
ature. Burrows fc Guol l| 1994 ) and lGuo fc Burrows! (|l997t > have reported 
on the ROSAT and ASCA X-ray observations of this remnant, finding a 
centrally peaked morphology which seems to be totally enclosed within 
the radio emission of this object. The latter, in turn, has a limb brightened 
double shell morphology, with a small diameter half-shell at NE (about 
30', a .k.a. the shell) and a larger incomplete shell at SW (a.k.a. the 
wing, iPineault et al. Hl985llPineault et al. 1 1 19871 . iLeahv fc Tian|[2005l and 
references therein). 

G166.0+4.3 was observed by XMM-Newton on 21 February 2003 (ID 
0145500101) and on 23 February 2003 (ID 0145500201). The first obser- 
vation was pointed at the western part of the remnant (the "wing"), for a 
total PN exposure after the screening procedure of 5.0 ks, while the sec- 
ond one was pointed to the small incomplete shell at the east. This latter 
pointing was severely contaminated by high energy proton flares and the 
net PN exposure time is only 1 ks, and we have not used it for spectral 
analysis. 

A mosaic map of this remnant in the band 0.3-5.0 keV, obtained us- 
ing the two XMM-Newton archive observations, is shown in Figure RJ 
which also incl ude the 1420 MHz radio map of the Canadian Galactic 
Plane Survey l|Tavlor et al J 12003). The X-ray morphology derived by 
previous observation s is confirmed, and th e bright region (named "west- 
ern bright knot" bv lBurrows fc Guolll994h is resolved for the first time, 
showing an elongation in the NW-SE direction. On the light of its ra- 
dio and X-ray morphology, this remnant may be considere d part of the 
sub-cl ass of mixed-morphology (MM) SNRs introduced by iRho fc Petrel 
(1 19981 ). The adopted explana tion of the peculiar radio morphology dates 
back to lPineault et ahl (1 19871 ) and involves the explosion in a moderately 
dense medium followed by the breaking out of the shock in a rarefied hot 
tunnel, followed in tur n by the interaction with again a denser medium. 
iGuo fc Burrowsl (|l997h have made the only measurements of metal abun- 
dances in the X-ray spectrum of G166.0+4.3, finding a substantial under- 
abundance of Mg, Si and Fe in the western bright knot and in the SE and 
NE of the wing. 

In order to study the abundances in the centrally peaked region of 
G166.0+4.3, we tried to characterize the thermal properties in the plasma 
of this remnant by selecting 4 interesting regions for spectral analysis, 
which ar e shown in Fig. HI Reg ions 1,2 and 3 are similar to region 1, 2 
and 3 of IGuo fc Burrowsl l|l997h . but they are considerably smaller and 
pointed to bright features, and therefore they are much less affected by 
background (also considering the smaller XMM-Newton PSF compared 
to the ASCA one ). In each region, we fitted the X-ray spectra with a 
iMewe et al.l (|l985l ) spectral model modified by interstellar absorption and 
letting the abundances of O, Ne, Mg, Si, S and Fe free to vary. In Fig. 
[5] we show, as an example, the spectrum of region 1. The results of the 
spectral analysis are presented in Fig. [S] and indicate low variation of 
temperature (up to 10%) and large variation of interstellar absorption 
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Figure 4: Top: XMM-Newton EPIC images of G166. 0+4.3 in the 0.3-5 keV 
background subtracted and vignetting corrected. Black thick contour is at 0.2 x 
10~ 3 cnt s _1 pix -1 and black thin contours correspond to 0.6, 0.8, 1.0, 1.2 and 
1.5 x 10~ 3 cnt s _1 pix -1 (1 pixel = 12" x 12"). Numbers mark the spectral 
extraction region used in the text (Region 1 is defined by the outermost thin 
surface brightness contour). The XMM-Newton EPIC field of view is marked 
in white (two observations combined ). Bottom 1420 MHz CGPS DRAO radio 



image of G166.0+4.3 (adapted from lLeahv fc Tianll2005f ) . X-ray contours of the 
top panel are shown in black. 9 
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Figure 5: XMM-Newton EPIC spectrum of the G166.0+4.3 Region 1 in Fig. H 
The best fit thermal model is shown as continuous line (with residual in lower 
panel). 

(up to 50%). The variation patterns see ms to be mo re in agreement 
with those found bv lBurrows fc Guol (|l994h than those of lGuo fc Burrows! 
(|l997l ). The regions seems to be characterized by an underabundance of 
O and Si and an overabundance of S, while Ne, Mg and Fe seems to be 
around solar values. O ur Mg, Si and Fe abundances are in agreement with 
iGuo fc Burrows! (|l997l ), who also found a value of < 1 for the abundances 
of these elements. A thermal model in Non-Equilibrium of Ionization do 
not improve the fits in the regions, giving a ionization time r > 10 12 
s cm -3 and abundances very similar to the one we obtained with the 
equilibrium model. 

We have investigated if there is a radial gradient of metal abundances 
variations inside the bright region 1 of the remnant (the "western bright 
knot"), in a manner similar to what we have done in Sect. 2 for IC443, 
but, given the low counting statistics, the results are not constraining. 

4 Discussion 

IC443 and G166.0+4.3 were listed in the lLazendic fc Slanel (|2006h work as 
standard abundances MMSNRs. However, the spatially resolved spectral 
analysis we have done on the bright X-ray central regions of IC443 has 
pointed out that this object belong to the subclass of metal rich mixed 
morphology remnant. In the case of G166.0+4.3 there is no compelling 
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Figure 6: Spectral fitting results of regions 1-4 in G166.0+4.3. See Fig. 
the position of the regions. Black arrows indicates a lower limit of 0. 



evidence for enhanced metal abundances apart from sulfur in the outer 
regions. The short exposure time of the XMM-Newton observation of the 
latter has prevented us to study possible trends of metal abundances in 
the core of this remnant . Our results suggest that it is important to review 
the list of MMSNRs of lLazendic fc Sland ((2006) to properly address the 
issue of the chemical abundances in all the known remnant of this kind. 
Other MMSNRs may be erroneously cataloged as standard abundances 
objects. 

It is not straightforward to understand if the high metal abundances 
seen in some of the MMSNRs is a general property, or if represent just an 
evolutionary phase of the life of these objects, or it depends on the initial 
remnant compositi on. From the theoreti cal point of view, the evaporat- 
ing cloud model of IWhite fc Lonel (|l99ll ) and the radiative SNR model 
I Cox et all l| 19991 ). traditionally used to explain the centrally peaked ther- 
mal X-ray morphologies, can be of limited help in understanding the sub- 
class of the metal rich MMSNRs, since they do not consider in detail the 
mixing between ejecta and shocked material inside the remnants. How- 
ever, these model ha ve been often used to interpret the X-ray emission o f 
some MMSNRs (e.g. iBocchino fc Bandierall2003l . lLazendic fc Slanell200fJ) . 
so it is worth to perform a detailed comparison with our results to under- 
line their limitations and to drive the exploration of new ones. 

We first consider the evaporating cloud model of lWhite fc Londi|l99lft . 
in which the X-ray central enhancement is due to the material ablated by 
shock heated clouds inside the remnant. The models depends on the tra- 
ditional E (the explosion energy), p (the preshock ISM density), t (the age 
of the remnant) plus two parameters describing the evaporating clouds, 
namely C (the ratio of mass in clouds to mass in the intercloud medium) 
and r (the ratio of the cloud evaporation time scale to the SNR age). 
With an appropriate choice of normalized quantities for the radial profiles 
(for what concern our study, we will consider the density normalized to 
the post-shock density, p/p s , the temperature normalized to the shell tem- 
perature T/T s , and the normalized emission measure EM{ivR 2 s )/EM(toi) 
see their Figure 4), they show that the dependence on E, po and t can 
be masked out in the presentation of the results. To ease the comparison 
with observations, the authors discuss the asymptotic behavior of their 
model when C and r — + oo, in which case the model only depend on the 
ratio C/t. In Fig. [7] we show the normalized temperature, density and 
emission measure for the models with C/t = 0, 1, 2, 3 and 4. To further 
ease the comparison with the data, we decided to normalize the temper- 
ature and density to their respective central valu es, instead of their shell 
values, like in the work of I White fc Longl (l991). In fact, in our obser- 
vations (and in general in all the MMSNR observations) the low X-ray 
surface brightness at the shell prevents an accurate spectral analysis, so 
it is particularly difficult to estimate the shell values. On the other hand, 
the central regions are by definition the brightest, and we may get more 
accurate results. One more point of concern in the comparison between 
models and our two targets is the value of the shell radius, used to nor- 
malize the abscissa of the profiles and the emission measure. IC443 and 
G166.0+4.3 have far from circular symmetric morphologies and offset X- 
ray peaks, so it is particular difficult to estimate a single shell radius. 
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Figure 7: Comparison between the IWhite fc Long! (|199ll ) evaporating cloud 
model and the observed profiles in IC443 and G166.0+4.3. Top: Normalized 
emission measure (see text for explanation). Middle: Density normalized to 
central value. Bottom: post-shock temperature normalized to central value. 
Black and red crosses are IC443 low and high temperature component data 
wit h uncertainties, respe ctively; blue curves are G166.0+4.3 data. Black curves 
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White fc Lonl (|l99ll ) profiles for models with C/t = (dashed), 1, 2, 3 and 



4. The distance from the center is normalized to the shell radius (assumed 10' 
for IC443 and - 6.5' for G166.0+4.3, see text for details). 
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For IC443, we have used a radius of 10' (~ 4.3 pc at 1.5 kpc distance), 
because it is the average distance of the radio shell from the X-ray peak 
in the East, North and West of the remnant. For G166.0+4.3, we have 
used ~ 6.5' (~ 13 pc at 4.5 kpc distance), which is the distance to the 
east radio wing. 

In Fig. [7] we also plot the derived values for the low and high tem- 
perature component of IC443 and for the single thermal component of 
G166.0+4.3. The profiles of G166.0+4.3 seems to be steeper then the 
profiles of IC443. In general, a reasonably good agreement can be found 
between the observed normalized emission measure and the model with 
C/t = 2 — 3 for both objects, even if the uncertainties about the adopted 
shell radius may affect the actual result. However, it is important to stress 
that the density and temperature profiles are not entirely consistent with 
the emission measure profiles. In case of IC443, the density profile would 
be more in agreement with C/t > 4, while the temperature profile with 
C/t < 1. It is worth noting that the inconsistency is present in both 
the thermal components of IC443, and that it also holds for the density 
of G166.0+4.3 which would indicate C/t > 4. Conversely, we could use 
the density profile to find a good fit to the model by varying the shell 
radius within a ±50% range, but the emission measure and temperature 
profile would be still inconsistent. For instance, we have verified that the 
emission measure would be much more centrally peaked than observed. 
Finally, we note that the central peak in both the observed density profiles 
see ms to exclude evapora ting cloud models with low values of r (cfr. Fig. 
1 in lWhite fc Lonslll99ll ). 

In summary, the comparison with the evaporating cloud model points 
toward high values of r, but cannot be considered very informative on the 
actual value of C and r, and indicates that the observed density profiles 
are usually steeper than what predicted by the model which best fit the 
emission measure profiles, a situation which cannot be explained in terms 
of an ejecta component in addition to the evaporating cloud component. 
In the latter case, in fact, we would expect the density profile be shallower 
than the model, because of the additional density provided by the ejecta. 
We therefore conclude that the X-ray morphology of these remnants is not 
compatible with the evaporating cloud model. The temperature profiles, 
though not entirely consistent with the model, seems to indicate in any 
case very limited variations, and therefore points toward efficient thermal 
conduction inside both the remnants. The failure of the cloud evapora- 
tion model and the high efficiency of the thermal conduction suggest an 
explanation in terms of the entr opy-mixed thermally conductive model 
suggested bv lShelton et al.l i|2004h to explain the MMSNR W44, possibly 
amended to take into account the peculiar environment in which these 
remnants are expanding. A detailed numerical model of the metal rich 
MMSNRs taking into account the mixing between ejecta and shocked ISM 
material, the thermal conduction and a realistic model for the environ- 
ment, is required to verify this scenario. 

The comparison with the radiativ e model of ICox et al.l (Il999f) can be 
done more straightforwardly following lLazendic fc Sland (|2006l ) assuming 
that the actual shock radius is the radius at which the remnant enters the 
radiative stage. This yields a lower limit of the ambient and central densi- 
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ties, following the relations of lCox et all l| 19991) . Using realistic radius val- 
ues, as the one we have used in the comparison with the cloudy ISM model, 
we get ambient density lower limits of the order of 10 cm -3 . While it is 
true that in both the remnants there is substantial evidence of expansion 
in dense environments, it is also true that it is very unli kely that the expan - 
sion occurred always in a medium so dense. For IC443, iTroia et all (120061 ) 
argued that the remnant evolved inside wi n d blo wn bubble, an hypothe- 
sis origin a lly su ggested by iBraun fc Stroml (QUI), while for G166.0+4.3, 
iGaenslerl (|l998l ) argued the remnant has expanded in a low density hot 
tunnel before to encounter recently a high density region. The hot tun- 
nel putative location is right where the X-ray emission peak is located. 
In both cases, therefore, it is very unlikely that a radiative model which 
assumes a long evolution inside a very dense medium is applicable, as it 
is whenever there is evidence of a massive progenito r star, which usually 
have strong pre-supernova winds. Ch en et all J2008) suggest that the X- 
ray emission of Kes 27 MMSNR may be explained in term of a shock 
reflected by a cavity wall, a mechanism which may be at work also in 
IC443 and G166.0+4.3, given the evidences on their environments. 

5 Summary and conclusion 

We report on XMM-Newton X-ray observation of the two supernova rem- 
nant IC443 and G166.0+4.3, with particular emphasis on the nature of 
their centrally peaked thermal X-ray emission above 1 keV. We confirm 
that the X-ray morphology of these remnants entitle them to be included 
in the class of mixed morphology SNRs. Contrary to what was previously 
known, we find that the chemical abundances of the bright central peak of 
IC443 are above solar, promoting this object in the recently pointed out 
subclass of metal rich MMSNRs. For G166. 0+4.3, there is no evidence 
of enhanced metal abundances, except maybe for sulfur in the outer re- 
gions, but the limited quality of the data prevented us to perform spatially 
resolved spectroscopy of the central peak. 

We derived profiles of metallicity, temperature, densit y and emission 
measu re, a nd we com p ared them with prediction of the IWhite fc Lond 
(|l99ll ) and ICox et all (1 19991 ) models, traditionally used to explain the 
MMSNR morphology in other objects. We have found that the observed 
profiles are inconsistent with the models, and we argue that a more de- 
tailed modeling including the mixing of ISM and ejecta mate rial is required 
to exp lain the observations. The entropy mixing model of IShelton et all 
(|2004l ) seems to be the most promising starting point for a more detailed 
modeling, though the effects of the complex environm ent in which the 
MMSNRs seems to be loca ted (ILazendic fc Slandl2006l ) and the projec- 
tion effects as discussed bv lPetrukl (|200lh should be taken into account. 
Temperatures, density, and chemical abundances profiles of MMSNRs are 
therefore the key observational constraints to compare with future mod- 
els. 

We thank D. Leahy for providing us with the 1420 MHz total in- 
tensity map of IC443 in electronic format. We also thank O. Petruk 
and S. Orlando for useful comments on this work. This work makes 
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